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@ System and method for a digital navigation satellite receiver. 

(57) There is disclosed a global positioning sys- 
tem receiver method and system which con- 
verts the analog signals to digital signals prior 
to perfonminng signal acquisition and GPS un- 
ique processing functions. The A/D converter 
(201) uses full null zone processing to increase 
anti-jamming capability and 4-level output to 
reduce signal processing. A. single analog and 
digital channel (302) is used for both the L1 and 
L2 channels as well as for all P(Y)-code and 
C/A-code thus allowing full multi channel track- 
ing with several channels (202, 203, 204) each 
tracking a separate signal while the remainder 
of the channels rove over all of the signals on 
the alternate L-band with programmable duty 
cycles. The system allows for Y-code substitu- 
tion for P-code by multiplexing a portion of the 
Y-code generated between multiple channels. 



O 
CM 
00 



Q. 

LU 







1-0 a 5= 
(yi L*j -cc 


5 ^ • — 1 


C2 ^ on 


is 


Q£ 3: -a: 




a- m 





m si ii §s 



<7i 

CD O 



ii 



Q£ — i S> 



Jouve, 18. rue Saint-Denis, 75001 PARIS 



EP 0 501 829 A1 



I 



TECHNICAL FIELD OF THE INVENTION 

This invention relates to navigation satellite receivers, and more particularly to a digital receiver having 
rapid signal acquisition. 

5 

BACKGROUND OF THE INVENTION 

The Navstar Global Positioning System (GPS) is used to determine exact geographic position (i.e., latitude, 
longitude and height above the earth) absolute velocity, as well as the exact time. The navigation device, 
10 receiver, must calculate the position velocity and the time by determining distance and relative velocity to a 
series of satellites. The velocity of the receiver is calculated from the doppler frequency shift of signals trans- 
mitted from space while the exact position of the receiver is calculated from the time shift of data due to the 
distance the signals must travel from the satellite. This distance is called range and the doppler shift yields range 
rate. 

15 A GPS receiver must receive signals generated from a satellite about 1 1 ,000 miles away. A GPS satellite 

transmits about 6 watts spread spectrum signal. The satellite and receiver, employ spread spectrum techniques 
to differentiate the signal from the noise. "Spread spectrum" means that as time passes, the frequency at which 
a signal is being transmitted will change or the instantaneous phase at which. the signal is being transmitted 
will change. Using conrelation techniques, the satellite receiver can match the spread spectrum signal coming 

20 from the satellite with an image of the signal that the receiver attempts to estimate. A precise match of the satel- 
lite's spread spectrum produces a potential signal processing gain of up to 53 decibels. The use of spread spec- 
trum techniques is essential to receive 6 watts of energy transmitted 1 1,000 miles away. At the antenna, the 
GPS signal is about 20dB below ambient cosmic noise. 

Using the GPS system, four transmitting satellites are required, to solve for the three spatial dimensions 

25 and time bias. The basic method of determining position is knowing the time difference from each of the satel- 
lites. The time difference for each satellite is the time required for a signal initiated at the satellite to be received 
by the user. Therefore, at least four satellites must be trackked. This is accomplished ideally by the use of at 
least four hardware receiver channels. 

Classically and historically, receiver channels have been large, power hungry and expensive. While this 

30 gives the best performance, it costs a lot of money. The lowest cost approach uses sequential tracking. Under 
sequential tracking, there is one receiver channel that sequences across multiple satellites-tracking each satel- 
lite for a predetemined period of time, and then tracking another satellite, etc. Sequential tracking requires the 
minimum hardware, but also has the lowest performance characteristics. 

There has been developed a multiplex technique which was essentially a compromise between using four 

35 receiver channels and using purely sequential tracking. This technique sequences very quickly across four 
satellites. This provides the dynamic capability of a continuous track receiver, but has a serious signal strength 
penalty under jamming conditions. The present state of the art requires a choice between either large and 
expensive receivers for high performance, or small and more affordable sequential tracking, with a severe com- 
promise in performance. 

40 Therefore, a need exists for a receiver that is small and affordable, yet with no compromise in performance. 

A further need exists in the art for a GPS receiver that achieves enhanced performance under diabolical 
conditions, such as jamming and fast acquisition. 

A still further need exists in the art for a receiver having simultaneous C/A and P or Y-code (P(Y)-code) 
capabality while reducing the parts count and enhancing signal acquisition time. 

45 

SUMMARY OF THE INVENTION 

The foregoing requirements have been met in this invention by a system and method which converts the 
analog signal to digital prior to processing the signal for detection purposes and which extends the prior use 
50 of channels on a chip. The receiver has rapid acquisition with high performance and anti-jamming capability 
with respect to the P(Y)-code as well as the C/A-code. 

Each GPS satellite generates two spread spectrum signals centered around separate frequencies. The L1 
channel is centered around 1575.42 megahertz, and has C/A-code, and P(Y)-code and data modulated on it. 
C/A-code has a 1.023 megahertz chipping rate with a bandwidth of about two megahertz, and P(Y)-code has 
55 a 10.23 megahertz chippinp rate with a bandwidth of about 20 megahertz. The L2 channel is centered around 
1227.6 megahertz and only has P(Y)-code and data modulated on it. Data provides information about the satel- 
lite's orbit and health and infomnation about other satellites in the GPS constellation. 

An existing digital receiver was designed using only C/A-code primarily because of semiconductor technol- 
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ogy that was available at that time. In order to design second generation or P(Y)-code receiver, advantage was 
taken of the higher signal processing speeds in silicon. One objective was to make the parts count lower while 
also achieving rapid signal acquisition. These objectives have been achieved by significantly reducing the pro- 
cessor loading requirements during signal search. In addition, the receiver was designed using modular con- 
5 cepts such that approximately 70 percent of the GPS receiver uses common circuitry for several receiver 
applications. 

In using the digital signal acquisition technique, advantage was taken of null zone automatic gain control 
in the analog to digital converter (A/D) to increase jamming immunity. The output from the A/D uses four level 
coding, thereby further reducing processing time. 
10 One technical advantage of our invention is that a GPS receiver is designed having a method and system 

for converting the analog signal to the digital equivalent prior to processing the signal for acquisition purposes. 

A further technical advantage of our receiver is that the processing time is reduced thereby allowing for 
faster signal acquisition during the signal search mode. 

A still further technical advantage of our invention is that the part count of the GPS receiver is reduced and 
15 the system is designed using modular blocks which are interchangeable with many different types of receivers 
and for receivers used for different applications. 

BRIEF DESCRIPTION OF THE DRAWINGS 

20 A more complete understanding of the present invention may be acquired by referring to the detailed des- 

cription and claims when considered in connection with the accompanying drawings in which like reference 
numbers indicate like features wherein: 

FIGURE 1 show a prior art navigation receiver; 
FIGURE 2 shows the navigation receiver of the present invention; 
25 FIGURE 3 shows a more detailed block diagram of one configuration of the present receiver; 

FIGURES 4a, 4b, 4c and 5 show details of null zone detector; 

FIGURE 6 shows one configuration of the receiver timing; ^ 
FIGURE 7 shows the chart of sensitivity by combining C/A and P(Y)-code measurements; 
FIGURE 8 shows a block diagram of two channel chips used in the receiver; 
30 FIGURE 9 shows a chart for discussion purposes relative to the receiver search and acquisition function; 

FIGURES 10a and 10b show details of the wide band noise meter; 

FIGURE 11 shows a block diagram of one embodiment of a discrete fourier transformer (DFT) system; 
FIGURE 12 shows complex vectors for the DFT circuit. 
FIGURE 13 shows a block diagram of the search processor; 
35 FIGURES 14, 15 and 16 show charts and data with respect to the discrete fourier transformer and the 

envelope detection; 

FIGURE 17 shows an algorithm for the basic Tong detection function; and 

FIGURE 18 shows a comparison chart of the various receivers and their relative performance. 

40 DETAILED DESCIPTION OF THE INVENTION 

FIGURE 1 shows prior art receiver 10 having antenna 12 which feeds into antenna preamp assembly 101. 
This is followed by L1/L2 frequency selector 102 and L-band down-conversion 103 which converts either the 
LI signal or the L2 signal to a second IF frequency. The local oscillators for the down-conversion are provided 
45 by reference oscillator 107 and frequency synthesizer 108. The second IF (output from conversion 103) goes 
into analog circuit 104 for code and carrier wipe off. This is where the code transmitted by the satellite is com- 
pared and matched with the receiver's estimate of the code, and the carrier is compared with the receiver's 
estimate of the carrier. 

Carrier frequency changes as a function of doppter velocity between the satellite and the user, and is the 
50 way the GPS measures relative velocity. Code generation 109 and canrier synthesizer 110 feed the code and 
carrier to the code/carrier wipe off circuitry and operate under control of computer 1 1 . 

After code and carrier wipe off, the signal is digitized in analog to digital (A/D) converter 1 05 and base band 
detected in module 106. The output of the base band predetection comprises the estimate of the error between 
the code and carrier that is transmitted from the satellite and the code and carrier prediction that is being made 
55 by the receiver. The en-or signal from the base band predetection to microprocessor 1 1 is used to generate 
new code and carrier states in modules 109 and 110. The key point is that significant signal processing is occur- 
ring in the analog domain, specifically code and carrier wipe off 104. 

For a more complete understanding of the typical receiver system reference is made to U.S. Patent 
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4,468.793 dated August 28, 1984, entitled "Global Position Systenn (GPS) Multiplexed Receiver," which refer- 
ence is hereby incorporated by reference herein- 
Turning now to FIGURE 2, there is shown receiver 20 which utilizes antenna 22 feeding into antenna 
preamp assembly 101 which in turn feeds into L1/L2 frequency selector 102. which in turn feeds into L-band 
5 down conversion 1 03, all operating under control of reference oscillator 1 07 and frequency synthesizer 1 08 in 
the same manner as similar numbered modules of receiver 1 0 of the prior art shown in FIGURE 1 . 

The output of L-band down conversion 103 is provided to analog to digital (A/D) AGC circuit 201 which is 
different from the A/D circuit in FIGURE 1. A/D circuit 201 employs a null zone Automatic Gain Control (AGC) 
technique which provides added immunity against jamming. The A/D converter provides a four level output sig- 
10 nal to signal processing circuits 202, 203 and 204. Note that circuits 202. 203 and 204 are replicas of the same 
chip and are detailed in concurrently filed copending patent application entitled "Method and System for a Dual 
Channel and Search Global Position System Signal Processor," which is hereby incorporated by reference 
herein. 

The outputs of circuits 202, 203 and 204 are shown as signals 205 going to processor computer 21 . They 

15 are similar to the outputs of FIGURE 1 going to processor 11 in that they represent the enror between the 
receiver's estimate of code and carrier phase and the actual code and carrier phase being received. Signal 
206 is the control from processor 21 back to circuits 202, 203, and 204, where the computer is directing the 
circuits to change their code and carrier estimates, as well as mode control of these chips. 

Modules 101, 102 and 103 contain circuitry common with that found in a good FM radio. At circuit element 

20 201 we have added null zone processing, which takes advantage of the spectra! characteristics of the GPS 
signal and the spectral characteristics of the most common GPS jamming signals, which are CW signals. The 
null zone circuitry is desensitized to process CW energy, thereby yielding up to a 9dB improvement in jamming 
immunity against CW jammers. 

Circuit elements 202, 203 and 204 channel signal processing with no signal differentiation occurring prior 

25 to these circuit elements. Each element uniquely and completely generates carrier estimates, code estimates, 
base band predetection estimates and contains all the correlators for signal processing, and provides auton- 
omous signal search capability for two satellites. Each circuit element also has the capability to simultaneously 
track both the C/A-code signal and the P(Y)-code signal within a single channel from a single satellite, effectively 
giving about a 4dB improvement in signal strength over tracking P(Y)-code only. 

30 FIGURE 3 shows a block diagram containing more detail of system 20 showing one implementation of the 

receiver. Note that modules 103 and 108 have been combined into a single block consisting of L-band hybrid 
circuit 302 and times 17 circuit 103. The exact arrangement of this circuit is dependent upon the arrangement 
of the various components, and this could be one circuit or multiple pieces separated and is driven by specific 
application requirements, 

35 It should be noted that reference oscillator 1 07 operates at an off beat frequency of 20.4608 megahertz or 

2 X 1 0.2304 megahertz. Conventional GPS reference frequency is 1 0.2300 megahertz. It is the frequency offset 
(in this case 400 HZ) that allows us to implement digital oscillators inside channels 202, 203 and 204. 

Processor 21 is an implementation of the computer using the Intel chip set such that element 304 is an 
Intel 80386 and element 305 is an Intel 80387 floating point unit. The computer also includes static RAM 306, 

40 EPROM 308 where the program code is stored, and interface and glue logic 307. Glue logic provided address 
decode functions, interrupt control functions, DMA functions, and miscellaneous housekeeping required for the 
computer. Element 31 provides the input output reference between the GPS core receiver and the rest of the 
worid. 

Element 31 has an I/O function 309, bus drivers 310 for parallel I/O and serial drivers 311 for a serial I/O. 
45 FIGURES 4-4C show aspects of the null zone A to D converter (201 , FIGURE 3). 

FIGURE 4A shows the functional implementation of the null zone technique with 4 bit outputs A, B, C and 

D. 

FIGURE 4B is a truth table for the outputs of FIGURE 4A, and FIGURE 4C shows the difference in spectral 
content between a gaussian wave form and a sinusoidal wave form. The sinusoidal wave fonm is what you see 
50 when a CW jammer is incoming along with the GPS signal. 

The advanced A/D conversion process shown in FIGURE 4A is a modified version of a scheme called null- 
zone detection and was adapted for use for two reasons: 

. Reduction in the number of digital data bits being processed, which substantially decrease the size and 
power consumption of the processor, making it more reliable and much less expensive; and 
55 . Additional resistance against constant envelope or continuous wave (CW) interference. CW jamming can 

be particulariy troublesome during signal acquisitions on the C/A-code. 

FIGURE 4A shows the different aspects of the A/D conversion process. Autommatic gain control (AGC) 
circuit 401 maintains the raw incoming analog signal at a constant level. Three threshold detectors 403-405 
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quantize the incoming voltage into one of four values. When the inconning wavefonm's probability distribution 
function is gaussian In shape, as would be the case for broadband noise, the four values are assigned weights 
of -16. -4, +4, +16. This Weighting approximates a linear detector. The automatic gain control (AGC) maintains 
the percentage of ±16 values as a constant. 

5 A separate detector 406 detenmines if the signal is. indeed, gaussian or if it has a strong CW component. 

The probability distribution function of these two cases is very different. If a strong CW component exists, then 
a waveform with equivalent power will spend a larger percentage of time near the high and low thresholds. The 
CW detector senses this condition, and the weights applied to the four values change to -16, -1, +1, +16. The 
±1 values deemphasize signal between the high and low thresholds where less GPS signal can be extracted. 

10 The waveform straddling the high and low thresholds is rich in GPS signal and maintains full weighting. 

FIGURE 5 shows the performance gain using this detector approach. At a CW jammer-to-signal ratio (J/S) 
of 50 dB, the detector using 4/16 weighting preserves 7 dB of additional signal, compared to a linear detector 
and 9 dB using the 1/16 weighting. The receiver uses this increase in signal power to operate despite much 
more CW interference. 

15 FIGURE 6 shows one common version of the receiver timing which demonstrates the flexibility of six - 

hardware channels by tracking combinations of LI and L2 simultaneously. This mode will have five channels 
tracking LI continuously and one channel roving on L2 across the five signals at least once every 30 seconds. 
By selecting the L1/L2 duty cycle in software, the receiver can easily be tailored for the unique conditions of 
future applications. 

20 FIGURE 7 shows the signal sensitivity gain by simultaneously tracking and vector summing the C/A and 

P(Y)-code signals from a single satellite. 

The C/A and P(Y)-code signals were combined inside the receiver to create 4.5 dB additional signal power 
when compared to P(Y)-code alone. The receiver can use this extra power against broadband jammers or 
beneath thick foliage. The satellite transmits the C/A-code with 3 dB more power than P(Y)-code to assist 

25 acquisition by receivers. But, before a conventional receiver can establish steady-state P(Y)-code track, it sus- 
pends tracking C/A-code. This receiver implementation removes that deficiency. 

FIGURE 7 illustrates the relative signal strengths of C/A, P(Y). and C/A plus P(Y)-code. Both C/A and P(Y)- > 
code signals are demodulated separately using the extra correlators in the PCOAC. Then software adds the 
two signals together and forms a composite vector sum. Because there is no correlation between the broadband ;^ 

30 noise of the two signals, the composite signal is 4.5 dB stronger than the P(Y)-code signal. 

Returning briefly to FIGURE 3, the receiver employs a fast multiplexing or switching technique between 
the L1 and L2 signals as shown. Note that both L1 and L2 enter hybrid circuit 302, but only one signal comes 
out. That is because the channels are capable of quickly switching between LI code state and carrier phase 
and L2 code state and carrier phase. This preserves a single signal path for LI and L2 having the effect that 

35 there is no interchannel bias between L1 and L2 and no interchannel bias between any of the channels. This 
technique also allows us to build a full function receiver with only a single RF signal path and a single A/D con- 
verter, thus eliminating at least 50 percent of the expensive RF circuitiy required for a GPS receiver. The RF 
switch is a conventional switch, but this technique is made possible because channels 202, 203 and 204 are 
capable of switching from a single code state and carrier phase to a different code state and carrier phase very 

40 quickly unlike conventional receiver architectures. 

FIGURE 8 is a block diagram of dual channel decoder 202, which is detailed in the above-identified copend- 
ing patent application entitled "Method and System for a Multi Channel and Search Global Position System 
Signal Processor." The details of operation of that system will not be repeated herein. 

FIGURE 9 represents pictorially the problem inherent with determining position from a series of satellites. 

45 There are two main requirements to track a satellite's transmitted signal. One is to detennine the code state 
(C/A-code or P(Y)-code) of the signal being received. The second is to determine the doppler effect (frequently 
offset) of the signal being received. The code changes are akin to two people talking next to each other. What- 
ever the speaker says, the hearer receives at the exact same time. However, in the speaker were to move away 
far enough, the speaker would output information, and that infomnation would arrive at the receiving person 

50 sometime later. Thus, if the speaker were reciting the alphabet, A, B, C. D in order, the receiver might be lis- 
tening to the letter ""A" while the sender is speaking the letter *'C". This is called the code being received. By 
knowing the delay, i.e., the number of bits different in the transmitted information, one could figure out the dis- 
tance, and thus the position of the receiver along a single straight line from the sender or satellite. This repetition 
from three other satellites yields precise position and time of the receiver. 

55 In addition, the relative velocity between the satellite and the receiver Is translated into a doppler shift in 

firequency of the transmitted canrier signal. Therefore, initial acquisition of a GPS signal from a single satellite 
is akin to a two dimensional problem where one dimension is the code or position state and the other dimension 
is the doppler or velocity state. The problem can then be described graphically in three dimensions (FIGURE 
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9). Samples of the incoming signal are integrated in time across the eight correlators within each channel. The 
results of the sampling are stored at intervals of JJ8 where T| is the total integration time. This is shown in the 
code/time plane. Then, via an eight point discrete fourier transform and adjacent doppler bin Interpolation, this 
is converted to 13 x 8 samples in the code/frequency domain. Thus in a single search time, the receiver can 
5 determine the existence or absence of signal in many more code and doppler states than are possible with a 
conventional receiver. 

FIGURE 9, for comparison, shows the search area 901' and 902' covered per channel and per integration 
time by two prior generation receivers 901 and 902. Note that the code/frequency area search for both Is a one 
dimensional search per channel and per integration time. This unique approach not only searches a larger area, 

10 but is a two dimensional search per channel per integration time. In other words, where a conventional receiver 
with eight correlators per channel would have looked through from two to eight possiblle locations of code and 
carrier, this receiver would have looked for 8 x 1 3 possible positions on the search plane/per channel. This can 
be likened to a child's battleship game, where a conventional receiver would have had up to eight chances to 
find the battleship while this receiver will have 8x13 chances to find the battleship each time it gets a turn. 

15 Returning briefly to FIGURE 8, noise meter 81 9 is tied to the ninth correlator in each channel. It has a high 

sampling rate tied to the DFT samples for search. It is programmable for C/A and P(Y)-code bandwidths and 
is capable of aiding search by setting the thresholds in the search processor for the Tong detectors. This is 
done without aiding or intervention or time line requirements from the computer. 

FIGURE 10A shows a typical noise meter, and FIGURE 10B shows the predetection integration sampling 

20 associated with predetection integrator 1005 for the P(Y)-code and C/A code for three different frequency plans. 
The noise meter measures the variance of integration samples. The noise meter consists of three major com- 
ponents: 

- a pre-detection integrator 1005 

- a power discriminator 1006 

25 - a post-detection integrator 1007. 

The pre-detection integration period determines the bandwidth over which noise power is measured and 
is normally set to roughly match the bandwidth of the code being received. Taking absolute values provides a 
discriminant for the variance of I and Q. Post-detection integration acts as an averaging process to improve 
the quality of the noise power estimate. The noise meter can measure the power of gaussian, CW or plus noise. 

30 The pre-detection integration period determines the noise power bandwidth. For nonnal operation, the 

bandwidth is set to match the bandwidth of the desired code (i.e., 10 MHz for P(Y)-code and 1 MHz for C/A- 
code). Pre-detection integration band width is equal to: 



pre-detection integration period 
35 For tracking operations, the noise meter post-detection integration period is the same as the sample inte- 

gration period. During search, the noise meter post-detection integration period is the same as the DFT period. 
Pre-detection integration samples are accumulated for the entire post-detection integration period. 
Variance of I and Q integration samples in a gaussian noise environment is calculated using the following 
equation: 

^ VAR,, = 0.0783321 x J/^^^.s ^ DATA„, [2] 

where VARfq is the variance of the I and Q integration samples, H is the pre-detection integration period, M is 
the number of reference clocks in the post-detection integration period, and DATAn^is the noise meter data. 
Note, there is a four reference clock dead time between pre-detection integration samples. 
45 The expected value of noise meter data in a gaussian noise environment is calculated using the following 

equation: 

DATAnm = 0.0249339 x ^^J^^ x SIG^^ [3] 
N + 4 

where DATAp^ is the expected value (mean) and SIGstd is the standard deviation of the incoming signal. Vari- 

50 ance of noise meter sample data is calculated using the following equation: 

VARnrr, = 0.363380 X X SIG^d2 [4] 

where VARn^ is the variance. 

For measuring gaussian noise, modulating P(Y)-code onto the incoming signal prior to pre-detection inte- 
55 gration has a negligible on noise meter data. However, if the noise has some CW content, modulating P(Y)-code 
onto the incoming signal prior to pre-detection integration reduces the effect CW noise has on final post-detec- 
tion values. The noise meter is also capable of measuring CW and/or pulse noise. 

FIGURE 8 shows the search engine works with the noise meter, the DFT, and the I and Q samples. First, 
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I and Q samples from the eight correlators are stored in I and Q sample buffering memory 821 . (I and Q means 
inphase and quadrature, which are the two phases which define a vector.) Once stored in memory, the DFT 
module 824 reads these samples out, applies the appropriate multiply vectors to each sample and adds the 
product to the accumulated value for each of the seven frequency bins. 
5 After a DFT cycle Is complete, the search processor 823 will go and read the samples and apply a Tong 

detect algorithm to them to detenmine If signal is present at a particular code state and carrier frequency. 

DFT module 824 Is shown in FIGURE 1 1 and uses an approximation of the standard eight-point DFT 
algorithm to process sample integration data into seven frequency bins based on a complex 2-bit vector set. 
As a pre-processor to search processor module 823 (FIGURE 8), DFT module 824 extends the frequency range 
10 of a search by up to a factor of seven, thus reducing search time in most applications. 

The approximate DFT algorithm differs from the standard eight-point DFT algorithm only by the vector set 
that is used. The approximate vector set is used to reduce the circuitry required for implementation. 

FIGURE 12 shows the vector set, [Vq, Vi, V2, V3, V4, V5, Ve, and V7], used for this approximate DFT 
algorithm. 

15 This algorithm can also expressed using the following matrix multiply equation: 
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where 2n represents the complex value In + jQn- Thus, a set of eight integration samples (Zq, Zi, Z2, Z3, Z4, Z5, 
Ze and Z7) are multiplied by a vector sequence producing frequency data for a given frequency bin. 

The output level associated with each doppler bin is affected by doppler attenuation. Attenuation increases 
45 as signal offset from the internal generated carrier center frequency increases. The result is a non-linear dis- 
tribution over the full range of the eight-point DFT output In addition to the seven frequency bins produced by 
the DFT module, the search processor module generates six frequency bins by interpolating between the seven 
bins generated by the DFT module. This provides a means of filling in the signal to noise ratio holes between 
frequency bins. 

50 Search processor 823 (FIGURE 8) controls the DFT circuitry. A discrete fourier transform is performed on 

the I and Q data sample integration set. Frequency data from the DFT consists of seven I and Q pairs (repre- 
senting seven doppler bins). The search processor interpolates between adjacent pairs of doppler bin data to 
obtain six additional doppler bins for a total of thirteen doppler bins for each of the correlators. A Tong detection 
algorithm is performed on this data to determine if a signal is present. If not, the search processor moves code 

56 to a new code state and continues the search there. The search continues until one of three conditions exist: 
a signal is found, the maximum search length (specified by the host processor) is reached, or an error condition 
occurs. Except for setup, the search processor operates independent of the resident computer. 

The seven frequency bins from the DFT are labelled Bo, B2, B4, Be. Bg, Bio and B12 from the most positive 
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to the most negative doppler, respectively. The search processor interpolates between these frequency bins 
providing bins B^, B3, 85. 87, 89, and B^^ These interpolated frequency bins Indicate signal power found at a 
center frequency which Is the arithmetic mean of the center frequencies of the even frequency bins from which 
It was interpolated. 

(l+jQ)n.1 = (l+jQ)n-(l+jQ)n.2 [6] 

where n = 0, 2, 4, 6, 8, and 10. 

Interpolation by way of the above equation means that odd doppler bins (interpolated values) will on aver- 
age have a magnitude 1 .414 times that of the magnitude of the even doppler bins. For this reason the search 
threshold used for Tong detection interpolated bins is 1.414 times the threshold for non-interpolated bins. 

Search processor 823 calculates envelopes for each complex (l+jQ) doppler bin value. Ideally, envelopes 
would be calculated as the square root of the sum of the squares; however, in orderto significantly reduce circuit 
complexity a Robertson approximation is used as follows: 

if [I] < [Q] then E = [Q] + [l]/2, [7] 
if [1] > [Q] then E = [I] + [Q]/2 

where [ ] indicates absolute value and E is the resultant envelope. 

Several factors of the configuration contribute to final envelope magnitude and standard deviation and 
channel C/No loss. 

FIGURE 14 shows gain of each DFT bin relative to the magnitude of the sample integration data. 

FIGURE 15 shows C/No loss for each doppler bin as a function of doppler. Theeffectof doppler attenuation 
is factored into the values shown on this FIGURE. Doppler attenuation is loss which occurs because the final 
signal IF does not exactly match the carrier generator digital oscillator frequency. 

FIGURE 16 shows the effect each stage of signal correlation has on C/No and envelope magnitude and 
standard deviation. These factors are useful for determining C/No loss, determining if envelopes will overflow, 
and for setting search threshold (Vti). 

For example, if carrier hard-limiting Is off, adjacent correlator presuming is on, the integration clock period 
is 250 microseconds, and the frequency of the incoming signal is 1,000 Hz off the carrier generator. (See 
copending patent entitled "Method and System for a Dual channel and Search Global Position System Signal 
Processor"). C/No loss is calculated as follows: 

-1.25 dB for adjacent correlator presuming 
-0,24 dB for envelope calculation 
+ -2,4 to -0.9 dB for doppler bin 10 (.25/250 /is = 1000 H2) 



-3.9 to -2,4 dB total C/No loss 
The maximum envelope magnitude can be calculated as follows: 

rms signal input voltage 

1.5 from the carrier generator 

1/512 from truncation 

14,448 from sample period (for reference clock = 57,792MHz) 
1.09 from envelope calculation 
X 16,1, the maximum signal gain of bin 10 



= 7^2,8 X rms signal input voltage 

This is valuable in detemnlning if envelope calculations will overflow for a given set of conditions. 
Noise variance will have a gain of: 
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SIG^^^, input noise voltage (1 sigma) 
1.5 from the carrier generator 
1/512 from truncation 

120.2 from sample period (for reference clock = 57.792 MHz) 
1.15 from adjacent correlator presuming 
1,12 from envelope calculation 
X 6.32 for even bins 



15 Thus, for this set of conditions, 2.87 is the ratio of the standard deviation of noise of the envelope to the 

standard deviation of noise at the input of the PCOAC. This factor is useful for calculating search threshold, 
Vti. 

The Vti and Vta thresholds are used In the Tong detection algorithm in detenmining whether or not signal 
is present as discussed previously. The host processor can either set Vti and W^z a* a fixed value or specify Vti 
20 and Vt2 to be calculated by multiplying the output of the noise meter by a specified gain (Gt). The reason there 
are two voltage thresholds is because the magitudes of odd doppler bins are greater than the magnitudes of 
even doppler bins by a factor of the square root of two, also as discussed previously. Vti is used for even doppler 
bins and Vta (Vti multiplied by the square root of two) for odd doppler bins. 

During search, the noise meter's post-detection integrations are continued until all DFT input samples are 
25 collected. The total integration period affects the magnitude of the final noise measurement. These equations 
were shown previously. 

Vti is always calculated at the completion of each DFT by multiplying noise measurement by a gain value 
(G+) specified by the host processor or it is provided by the host processor directly. 
The search processor calculates Vt2 by multiplying Vti by the square root of two. 

30 A diagram of the Tong detection algorithm is shown in FIGURE 1 7. The Tong detection algorithm is desig- 

ned to decide when signal is present. The Tong detector keeps a Tong count for each of the 104 (13 doppler 
bins times 8 correlators) search bins for which it is trying to determine if signal is present. Initially, the Tong 
count (K) is set to one 1712. The search processor 1703 compares each envelope 1702 to its appropriate 
threshold (Vti or Vt2). If the envelope is greater than or equal to the threshold, the Tong count for that search 

35 bin is Incremented 1 704. If the envelope is less than the threshold, the Tong count is decremented 1 705. When 
a Tong count for a search bin reaches zero 1707, the signal at that search bin is declared "not found" and that 
search bin is no longer active. (Tong detection is no longer performed on it). If all 1 04 Tong counts reach zero, 
the signal is declared "not found" for all search bins at the current code state. When this occurs, the code is 
retarded 1710 to a new code state. The Tong detector is reset 1712 (all Tong counts are set to 1) and search 

40 continues. Signal is declared "found" 1708 if a Tong count is incremented to the value "A" 1706 set by the host 
processor. 

For some applications it may be desirable for the search processor to ignore the envelopes from certain 
correlators or for certain frequency bins. The search processor provides the host processor the means with 
which to mask the unwanted correlators or frequency bins. 

45 A search that is in progress can be temporarily suspended. This might be done in the event that the host 

processor has determined the input signal to be no good, such as when a pulse jammer is detected. Search 
suspend is tied to DFT cycles. Tong detection is not perfomned on DFT data which was taken during search 
suspend. Once suspend is de-asserted, Tong detection resumes on data from the next DFT cycle. 

Each channel has a (hardware) suspend search input. The host processor enables and disables the sus- 

50 pend input for a channel by setting a suspend input enable bit. In addition, the host processor can directly sus- 
pend search. 

The search processor implements a mush counter which counts the number of DFT periods the search 
processor has spent on a particular code state. The purpose of the mush counter Is to prevent the search pro- 
cessor from getting hung up at a particular code state because of a bad threshold, a jamming signal or high 
55 side-lobe. If the mush counter determines that the number of DFT samples spent on a particular code state is 
equal to a maximum mush count specified by the host processor, then the mush counter function can either 
halt search or force the search processor to go on to the next code state, depending on a bit set by the host 
processor. 
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The search processor also incorporates a hit counter which can be used by the host processor to aid in 
determining the adequacy of a threshold setting. It can also play a role in designs for non-gaussian environ- 
ments. The hit counter only provides information to the host processor; it does not affect the operation of the 
search processor. The hit counter counts the number of hits for a specified number of DFT samples. It is incre- 
5 mented each time in active Tong count is incremented. Once the hit count is complete, it is transferred to a 
latch to be read by the host processor. The hit counter is cleared after the transfer to keep count of hits for the 
next set of DFT samples. A large hit count could be indicative of a low threshold. 

FIGURE 18 summarizes the perfomnance of four different receiver architectures under various conditions 
of signal aiding and jamming levels. The metric used is time to first fix or the time from the point at which the 
10 receiver is turned-on until it is navigating to specified accuracy. In all cases the unique receiver architecture 
which is P(Y) described herein produces a system that gives the best performance independent of jamming 
level and independent of aiding quality. 

It should be noted that the GPS Satellite System has in it an algorithm which converts the P(Y)-code to 
Y-Code for special applications. This algorithm is not widely known. Therefore only those receivers that have 
15 a special access to the algorithm can utilize Y-code. 

In operation a special code called the CV code can be sent to the satellite, and the satellite using its special 
algorithm, will then convert the P(Y)-code signals to a Y-code signal. Receivers which have been outfitted with 
a special circuit which is available from the U.S. government, also in response to a corresponding CV code, 
will convert to a Y-code receiver. In this context CV code stands for crypto variable. In prior receivers, there has 
20 been one AOC circuit for each channel. In the instant disclosure there is one AOC alternative circuit embedded 
in each dual channel chip. Therefore, the AOC alternative each services two channels. Accordingly, once the 
CV code has been received or inserted into the receiver, the receiver becomes a six channel Y-code receiver 
with only three AOC altemative circuits embedded therein. The AOC alternative, which is used in place of the 
standard AOC circuit available from the govemment, takes advantage of the fact that the circuitry is shared 
25 between two channels and that by being embedded in the same chip with the search processor it can support 
the high code chipping rates of 28 megahertz and also the stopping of the code associated with search code 
state advances and retards. We do this by multiplexing a portion of the circuit. 

Although this description describes the invention with reference to the above specified embodiments, the 
claims and not this description limit the scope of the invention. Various modifications of the disclosed embodi- 
30 ment, as well as alternative embodiments of the invention, will become apparent to persons skilled in the art 
upon reference to the above description. Therefore, the appended claims will cover such modifications that fall 
within the true scope of the invention. 

35 Claims 

1. A receiver for determining a global positioning parameter from analog signal information received from a 
plurality of navigational satellites, said signal information having data shifts dependent upon distance and 
frequency shifts dependent upon the relative velocity between the satellite's transmission position and the 
40 geographical position of said receiver, and wherein said analog signal is significantly weaker than ambient 

atmospheric noise, said receiver comprising: a digital signal processing circuit for acquiring valid signals 
in the presence of said noise, for estimating the correct position and for correcting said estimate based 
upon data exchanged between said acquisition circuit and said processor. 

45 2. A receiver for determining a global positioning parameter from analog signal information received from a 
plurality of navigational satellites, said signal information having data shifts dependent upon distance and 
frequency shifts dependent upon the relative velocity between the satellite's transmission position and the 
geographical position of said receiver, and wherein said analog signal is significantly weaker than ambient 
atmospheric noise, said receiver comprising: 

so a first stage for accepting said analog signals from a plurality of said satellites and for modulating said 

received analog signals into equivalent IF analog signals; 

an analog to digital converter for converting said modulated IF signals into digital signals prior to further 
processing of said IF signals; 
a computer, and 

55 a digital signal processing circuit for acquiring valid signals in the presence of said noise, for estimating 

the correct position and for correcting said estimate based upon data exchanged between said tracking 
circuit and said computer. 
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The receiver set forth in claim 2 wherein said converter includes a full null zone detector. 

The receiver set forth in claim 2 wherein said digital signal processing circuit includes a noise meter for 
aiding In said signal acquisition or aiding In signal tracking. 

The receiver set forth in claim 2 wherein said analog satellite signal information is coded in P(Y)-code and 
in C/A code, and wherein said digital signal processing circuit includes: 
circuitry operable to track both said information codes simultaneously on a single channels. 

The receiver set forth in claim 2 wherein said receiver Includes an offbeat oscillator frequency which allows 
for digital oscillator to resolve carrier and P{Y)-code clock phase in order to process said P(Y)-code. 

7. The receiver set forth in claim 5 wherein said digital signal processing circuit includes: 

a plurality of signal correlators for aiding In said signal acquisition of P(Y)-code and C/A code. 

15 8, The receiver set forth in claim 1 further comprising: 

circuitry for tracking multiple channel is simultaneously, said circuitry Including: 
circuitry for assigning a different signal to a plurality of said channels; and 
circuitry for time sharing a plurality of signals on the remainder of said channels. . 

20 9. The receiver set forth in claim 2 wherein said analog satellite signal information Is coded in P(Y)-code and 
in C/A code and wherein said receiver further Includes: 
circuitry for receipt of satellite data indicating status of anti-spoofing ; and 

circuitry operative in response to receipt of said status for converting said receiver from P-code operation 
to Y-code operation. 

25 

10. The receiver set forth in claim 9 wherein said last mentioned converting circuitry Is shared between a 
plurality of channels. 
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11. A method for determining a global positioning parameter from analog signal Information received from a 
plurality of navigational satellites, said signal information having data shifts and frequency shifts dependent 
upon the distance and relative velocity between the satellite's transmission position and the geographical 
position of said receiver, and wherein said analog signal is significantly weaker than ambient atmospheric 
noise, said method comprising the steps of: 

accepting said analog signals from a plurality of said satellites and modulating 
said receiver analog signals Into equivalent IF analog signals; 

converting said analog IF signals to digital signals prior to further processing of said IF signals; 
processing said digital signals for acquiring valid signals In the presence of said noise; 
estimating the correct position; and 

continually updating said estimate based upon data exchanged between said estimating function and said 
signal processing function. 

12. The method set forth In claim 1 1 wherein converting includes the step of: 
Implementing full null zone detection. 

13. The method set forth in claim 1 1 wherein said analog satellite information Is coded In P(Y)-code and in 
C/A code, and wherein said signal processing function includes: 

tracking both said Information codes simultaneously on a single channel. 

14. The method set forth in claim 1 1 wherein said estimating function includes: 
vector summing said acquired P(Y) and C/A codes. 

50 

15. The method set forth in claim 1 1 wherein said signal processing function includes 

the use of an offbeat oscillator frequency which allows for phase resolution of digital oscillators to process 
said P(Y)-code. 

55 16. The method set forth in claim 1 1 further comprising the function of: 
tracking multiple channels simultaneously, said tracking including: 

assigning a different signal to each of several said channels; and time sharing a plurality of signals on said 
remainning channel. 
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17. A system for acquiring signals from at least one global positioning system satellite, said system comprising: 
and analog receiver section for receiving said signals together with atmospheric noise and for developing 
an IF signal containing a spread spectrum of signals; 

circuitry for converting said IF signal to a digital equivalent; and 
circuitry for locking on to valid signals, said locking circuitry including: 
a digital signal processor; and 

a host computer for control and processing of data contained in said signal; said 

data containing information for allowing estimates to be made with respect to the frequency and data shifts 
of said signal. 

18. The system set forth in claim 17 further comprising: 

circuitry for selectively adapting any channel on any of said signal processing circuits to process either or 
both P(Y) and C/A signals from said satellite. 

19. The system set forth in claim 17 wherein each channel of said multi-channel signal processing circuit is 
operable for receiving both P(Y)-code and C/A code simultaneously, said system further including: 
circuitry for vector summing said P(Y) and C/A code to aid in tracking said signal. 

20. The method for acquiring signals from at least one satellite of a global positioning system, said method 
comprising the steps of: 

receiving said signals in analog form together with atmospheric noise and for developing an IF signal con- 
taining a spread spectrum of signals; 
converting said if signals to a digital equivalent; and 

locking onto valid signals, said locking step including the reiterative step of: 

estimating received signals with respect to frequency and code state of said signals with a processor, said 
estimates being compared with actual signals for allowing new estimates to be calculated so as to perform 
said locking. 

21. The method set forth in claim 20 wherein said converting step includes the step of: 
performing null zone detection. 

22. The method set forth in claim 20 wherein said method includes the step of: 
using an off beat oscillator frequency for resolving signal phase. 

23. The method set forth in claim 20 wherein said locking step includes the use of a plurality of dual channel 
signal processing circuits each operable to perform said comparison. 

24. The method set forth in claim 20 further comprising the step of: 

selectively adapting any channel on any of said processing circuits to receive P(Y) or C/A signals simul- 
taneously from said satellite. 

25. The method set forth in claim 20 wherein each channel said dual channel signal processing is operable 
for receiving both P(Y)-code and C/A code simultaneously, said method further including the step of: 
vector summing of said P{Y) and C/A code to aid in tracking said signal. 

26. A receiver as claimed in claim 1 or claim 2 and wherein said parameter is position, velocity or time. 

27. A method as claimed in claim 1 1 wherein said parameter is position, velocity or time. 
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FIG. 6 
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FIG. 1 6 
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